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INTRODUCTION  fracture  toughness  with  increasing  temperature, 

i.e.,  a  brittle-to-ductile  transition  [4].  The  theory  is 
This  work  investigates  the  dynamic  response  implemented  in  a  dynamic  Lagrangian  finite 

of  a  two-phase  tungsten  alloy  of  interest  to  the  element  (FE)  code  [5],  and  2D  FE  meshes  are 

defense  industry  for  use  in  Kinetic  Energy  (KE)  constructed  of  the  initial  microstructures  of  a  KE 

penetration  applications  [1].  The  material  consists  rod  sectioned  either  transversely  or  parallel  to  the 

of  pure  tungsten  (W)  grains  of  body-centered-cubic  rod  axis  (i.e.,  the  extrusion  direction).  Simulations 

structure  embedded  in  a  face-centered-cubic  binder  of  high  rate  impact,  wave  propagation,  and  spall 

phase  of  W,  nickel  (Ni),  and  iron  (Fe).  In  what  are  conducted,  reminiscent  of  conditions  present  in 

follows,  we  develop  crystal  plasticity-based  plate  impact  experiments  [6,  7].  We  show  how 

constitutive  models  for  each  phase  that  account  for  statistical  data  acquired  from  simulations  depict 

the  effects  of  initial  lattice  orientation  (i.e.,  texture)  effects  of  the  microstructure  and  model  parameters, 

on  the  anisotropic  stress-strain  response  [2,  3]. 

Additionally,  cohesive  zone  models  are  formulated 

that  account  for  separation  at  grain  and  phase  THEORY 

boundaries  as  a  result  of  the  resolved  stress 

components  acting  on  these  boundaries,  and  Only  essential  aspects  of  the  constitutive 

reflecting  the  experimentally  observed  increase  in  theory  are  given  here;  the  interested  reader  may 
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find  more  details  in  the  references  [2,3].  The 
kinematics  of  thermo-elastoplasticity  within  each 
phase  of  the  bulk  material  is  described  by  the 
following  decomposition  of  the  deformation 
gradient  f: 


f  =  dx/dX  =  fef0fp ,  (1) 


where  fe  accounts  for  rotation  of  the  lattice  and 
stretching  due  to  applied  stresses,  fe  reflects 
stress-free  isotropic  thermal  expansion/contraction, 
and  f p  is  the  plastic  stretch  and  rotation  due  to 
inelastic  mechanisms,  namely  dislocation  glide. 
Spatial  and  referential  coordinates  are  denoted  by  x 
and  X,  respectively.  The  free  energy  per  unit  mass 
\j/  of  each  crystal  is  decomposed  as 

!¥w  =\-K0{d)32 -\k^  + 

/  (2) 

»{oyi*'’+-Kn{o)?+y(e), 


where  p  is  the  spatial  mass  density,  c  is  the 
specific  heat  capacity,  r(a)  is  the  resolved  Cauchy 
stress  on  slip  system  a ,  ^  is  the  plastic  slip  rate, 

and  k  is  the  thermal  conductivity.  The 
viscoplastic  flow  rule  is  given  by 


*) 


a(«) 

J 


sgn 


(4) 


with  the  resolved  Kirchhoff  stress,  g(a)  the 
evolving  temperature-dependent  slip  resistance, 
and  and  m  material  parameters. 

Interface  separation  is  modeled  by  the  cohesive 
zone  approach,  whereby  crack  initiation  occurs 
when  the  resolved  local  stress  component(s) 
exceed  a  user-prescribed  value.  In  the  present 
work,  fractures  occur  when  one  of  the  following 
criteria  is  achieved  locally: 

S=S0+Sx(6-90),  T  =  T0+Tl(d~e0),  (5) 


where  ft/c  is  the  mass  density  in  the  intermediate 

configuration,  6  is  the  temperature,  &  and  ee'  are 
respectively  the  scalar  elastic  volumetric  strain  and 
the  elastic  deviatoric  strain  tensor  in  the 
intermediate  frame,  %  is  an  internal  variable 
associated  with  the  internal  energy  contribution 
from  the  microscopic  stress  fields  of  lattice  defects, 
and  y  accounts  for  the  specific  heat.  The  nominal 

bulk  modulus  and  shear  modulus  are  given  by  K0 
and  p  ,  respectively,  a:  is  a  scalar  parameter,  and  a 
Kx  leads  to  a  volume-  (i.e.  pressure-)  dependent 
effective  bulk  modulus.  Note  that  W  is  virtually 
elastically  isotropic  at  room  temperature  and  low 
pressure.  The  local  energy  balance  is  written  as 


temperature 

change 


plastic 

dissipation 


energy  of 
lattice  defects 


(3) 


+HFS’ 


thermoelastic  coupling 


heat 

conduction 


where  s  is  the  resolved  nominal  normal  stress  at 
the  interface,  f  is  the  resolved  nominal  shear 
stress  at  the  interface,  and  the  remaining  symbols 
entering  Eq.  (5)  are  material  parameters.  Upon 
damage  initiation,  the  response  of  the  degraded 
material  at  interfaces  is  dictated  by  the  following 
irreversible  traction-displacement  relations: 


.?  =  (.?o +■?,(#- 00  )) 

i  =  (5o+ii(0-0o)) 


f  S  ^ 
1 — i 

v 


(load) 


(6) 


(unload) 


V  "c  J  "nl 


:  =  (fo  +  f,(0-6>o))x 


i-M 

v  SeJ 


\  M 

v  0,  y 


sgn  (A),  (load) 


:  =  (f0+f1(6'-6>0))x 


(7) 


V  0c  , 


'c  J 


K 


(unload) 
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with  8n  the  normal  crack  opening  displacement, 
8nX  the  maximum  previous  value  of  Sn,  St  the 
tangential  opening  displacement,  Sn  the  maximum 
previous  value  of  8t ,  and  8c  a  material  parameter 
defining  the  separation  distance  beyond  which  the 
damaged  zone  behaves  as  two  free  surfaces. 

NUMERICAL  IMPLEMENTATION 


The  constitutive  models  for  the  pure  W  and 
binder  phases,  as  well  as  those  for  cohesive 
interfaces,  were  implemented  within  the  EPIC 
dynamic  finite  element  program  [5].  For  details 
regarding  the  FE  implementation,  please  see  [3]. 

Unstructured  finite  element  meshes  conforming 
to  realistic  microstructures  were  constructed  from 
optical  images  of  tungsten  alloy  KE  rods,  as 
depicted  in  Fig.  1  and  Fig.  2.  Initial  lattice 
orientations  of  the  grains  were  designated  to  match 
experimentally-obtained  data  from  an  electron 
backscattering  diffraction  technique.  The  domain 
size  was  L  =  150  jum ,  and  on  the  order  of  20000 

plane  strain  triangular  elements  were  used,  with 
more  refined  areas  of  the  mesh  concentrated  in  the 
vicinity  of  grain  and  phase  boundaries.  The  OOF 
meshing  software  [8]  was  used  to  this  end.  Four- 
node  cohesive  elements  were  inserted  along  these 
boundaries,  with  coincident  nodes  constrained  to 
remain  rigidly  bonded  in  the  simulations  until 
either  of  the  criteria  (5)  was  attained. 


X2 


X1 


Figure  1.  Transverse  section  of  KE  rod  and  FE  mesh  of 
corresponding  equiaxed  granular  microstracture. 


Figure  2.  Longitudinal  section  of  KE  rod  and  FE  mesh 
of  corresponding  elongated  granular  microstructure. 


RESULTS  AND  DISCUSSION 

Here  we  report  on  impact,  wave  propagation, 
and  spall  behavior  of  the  tungsten  alloy.  For  these 
simulations,  the  following  boundary  conditions 
were  applied: 

•  Along  X1  =  0 ,  a  short  velocity  impulse  of 
duration  10  ns  was  applied  in  the  positive  1- 
direction,  corresponding  to  an  impact  velocity 
in  a  plate  impact  experiment  of  250  m/s. 
Subsequently,  the  boundary  X1  =  0  was  held 
rigid. 

•  Along  X1  =  L  ,  null  traction  conditions  were 
specified. 

•  Along  X2  =  0  and  X2  =  L ,  either  non- 
reflecting  conditions  [5]  or  uniaxial  strain 
conditions  (i.e.  null  velocity  in  the  2-direction) 
were  imposed  in  various  simulations. 

•  Null  heat  flux  conditions  were  specified  along 
all  external  surfaces. 

The  above  boundary  conditions  create  a  rightward- 
moving  compressive  wave  that  travels  through  the 
grid,  reflects  off  the  free  surface  at  X1  =  L  as  a 
tensile  wave,  then  potentially  leading  to  spall 
fracture  within  the  sample. 

Deformed  FE  meshes  for  an  equiaxed 
microstructure  in  Fig.  3,  at  an  elapsed  simulation 
time  of  100  ns,  illustrate  the  spall  separation  at  the 
material  interfaces.  On  the  left  image  in  Fig.  3,  we 
used  Sl=fl=  0  in  Eqs.  (5)-(7),  meaning  that  the 
brittle-ductile  transition  was  not  addressed.  On  the 
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right  image  of  Fig.  3,  nonzero  values  were  used  for 
sl  and  f, .  Different  fracture  patterns  emerge,  with 
fracture  suppressed  in  the  right  figure  relative  to 
the  left  due  to  cohesive  strengthening  following  the 
temperature  rise  that  occurs  commensurate  with 
wave  propagation  and  plastic  dissipation. 


Figure  3.  Deformed  FE  meshes  at  t  =  100  ns. 
Temperature-independent  cohesive  law  used  on  left 
simulation,  temperature-dependent  law  used  on  right. 


Figure  4.  Standard  deviation  in  free  surface  velocity  for 
selected  simulations  featuring  various  microstructures 
and  cohesive  parameters. 

Figure  4  shows  the  time  history  of  the  standard 
deviation  in  the  free  surface  velocities  recorded  at 

X1  =  L  .  Prior  to  t  &  25  ns  ,  velocities  are  zero  as 
the  compressive  wave  has  not  yet  reached  the 
surface.  We  see  that  the  velocity  statistics  are 
influenced  by  grain  shape  (i.e.,  equiaxed  versus 
elongated)  and  choice  of  cohesive  parameters.  In 
the  extreme  case  of  infinite  cohesive  strength  in 
which  spall  fractures  are  suppressed,  the  standard 
deviation  values  are  generally  much  smaller  than 
the  other  cases  in  which  spallation  occurs. 


CONCLUSIONS 

A  numerical  model  of  plastic  deformation  and 
spallation  in  multiphase  metallic  crystals  has  been 
exercised,  capturing  finite  anisotropic  elasto- 
plasticity,  thermodynamics  of  energy  storage, 
dissipation,  and  heat  conduction,  as  well  as 
interfacial  fracture.  The  results  presented  here 
indicate  relationships  between  microstructural 
features,  fracture  patterns,  and  free  surface  velocity 
statistics.  Such  information  should  provide 
increased  understanding  of  experimentally- 
obtained  results  from  plate  impact  tests  [6,  7]. 

ACKNOWLEDGEMENTS 

Support  from  the  Weapons  and  Materials 
Research  Directorate  of  the  U.S.  Army  Research 
Laboratory  is  acknowledged. 


REFERENCES 

1.  Magness,  L.S.,  “High  strain  rate  deformation 
behaviors  of  kinetic  energy  penetrator  materials 
during  ballistic  impact”,  Mech.  Mater.  17,  147-154, 
1994. 

2.  Clayton,  J.D.,  “Dynamic  plasticity  and  fracture  in 
high  density  poly  crystals:  constitutive  modeling  and 
numerical  simulation”,  J.  Mech.  Phys.  Solids  53, 
261-301,2005. 

3.  Clayton,  J.D.,  “Modeling  dynamic  plasticity  and 
spall  fracture  in  high  density  polycrystalline  alloys”, 
Int.  J.  Solids  Structures  42,  4613-4640,  2005. 

4.  O’Donnell,  R.G.  and  Woodward,  R.L.,  “Influence 
of  temperature  on  the  fracture  of  a  W-Ni-Fe  alloy”, 
J.  Mater.  Sci.  35,  4319-4324,  2000. 

5.  Johnson,  G.R.,  Stryk,  R.A.,  Holmquist,  T.J.,  and 
Beissel,  S.R.,  “Numerical  algorithms  in  a 
Lagrangian  hydrocode”  ,  Wright  Laboratory, 
U.S. A.,  WL-TR- 1997-7039,  1997. 

6.  Vogler,  T.J.,  Reinhart,  W.D.,  and  Chhabildas,  L.C., 
“Dynamic  behavior  of  boron  carbide”,  J.  Appl. 
Phys.  95,4173-4183,2004. 

7.  Vogler,  T.J.  and  Clayton,  J.D.,  “Heterogeneous 
spall  failure  of  a  tungsten  alloy”,  Proc.  Soc.  for 
Experimental  Mechanics  Annual  Conference  and 
Exposition,  Portland,  OR,  USA,  June  7-9,  2005. 

8.  Langer,  S.,  Reid,  A.,  Carter,  C.,  Fuller,  E.,  and 
Roosen,  A.,  PPM200F  vl.l ,  freeware  website 
http://www.ctcms.nist.gov/,  2003. 


314 


NO.  OF 

COPIES  ORGANIZATION 


1  DEFENSE  TECHNICAL 
(PDF  INFORMATION  CTR 
ONLY)  DTICOCA 

8725  JOHN  J  KINGMAN  RD 
STE  0944 

FORT  BEL  VOIR  VA  22060-6218 

1  US  ARMY  RSRCH  DEV  & 
ENGRG  CMD 
SYSTEMS  OF  SYSTEMS 
INTEGRATION 
AMSRD  SS  T 
6000  6TH  ST  STE  100 
FORT  BEL  VOIR  VA  22060-5608 

1  DIRECTOR 

US  ARMY  RESEARCH  LAB 
IMNE  ALC  IMS 
2800  POWDER  MILL  RD 
ADELPHI  MD  20783-1 197 

3  DIRECTOR 

US  ARMY  RESEARCH  LAB 
AMSRD  ARL  Cl  OK  TL 
2800  POWDER  MILL  RD 
ADELPHI  MD  20783-1 197 


ABERDEEN  PROVING  GROUND 
1  DIR  USARL 

AMSRD  ARL  Cl  OK  TP  (BLDG  4600) 


NO.  OF 

COPIES  ORGANIZATION 


ABERDEEN  PROVING  GROUND 


32  DIR  USARL 

AMSRD  ARL  Cl  HC 
P CHUNG 
AMSRD  ARL  WM 
J  MCCAULEY 
T  WRIGHT 

AMSRD  ARL  WM  MA 
W  NOTHWANG 
AMSRD  ARL  WM  TA 
S  SCHOENFELD 
AMSRD  ARL  WM  TC 
M  FERMEN  COKER 
R  COATES 
AMSRD  ARL  WM  TD 
S  BILYK 
T  BJERKE 
D  CASEM 
J  CLAYTON  (5  CPS) 
T  CLINE 
D  DANDEKAR 
W  EDMANSON 
M  GREENFIELD 
C  GUNNARSSON 
Y  HUANG 
K  IYER 
R  KRAFT 
BLOVE 
S  MCNEILL 
H  MEYER 
RMUDD 
M  RAFTENBERG 
E  RAPACKI 
M  SCHEIDLER 
S  SEGLETES 
T  WEERASOORIYA 


